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4th-0RDER  TENSOR  INVARIANTS  AND  GEOMETRIC  REPRESENTATION 

Edward  M,  Wu* 

ABSTRACT 

First  and  second  order  invariants  of  4th-order  tensor  are 
derived.  Geometric  representations  analocpus  to  Mohr's  circle 
are  presented  to  aid  the  visualization  and  operations  of  4th 
order  tensors  such  as  transformation,  determination  of  principal 
direction  and  optimization  procedures.  Possible  applications  of 
the  invariants  are  suggested. 


♦Assistant  Professor,  Materials  Research  Laboratory,  Mechanical 
and  Aerospace  Engineering,  Washington  University,  St.  Louis, 


Missouri . 


Introduction 


Many  physical  properties  of  composites  and  crystals  are 
functions  of  material  orientations.  Such  properties  can  be 
characterized  by  the  appropriate  constitutive  functions.  In 

order  to  insure  the  invariancy  with  respect  to  coordinate 
transformations,  these  functions  are  customarily  expressed  in 

tensorial  forms.  Additional  groundwork  is  required  in  practical 
engineering  application  of  such  constitutive  functions.  Material 
constants  must  be  measured  for  design  computations,  and  analytical 
techniques  must  be  explored  for  operational  efficiency.  Both 
experimental  techniques  and  operational  procedures  are  well 

established  for  physical  properties  which  are  scalars,  vectors 
and  second  order  tensors.  For  example,  the  measurement  of  strain 
tensor,  as  well  as  its  analysis,  visualization  and  transformation 
by  the  Mohr's  Circle  are  well  known.  Howaver,  comparable 
operations  for  4th- order  tensors  have  not  been  fully  explored.  In 
the  advent  of  composites  as  anisotropic  engineering  materials  , 
the  analysis  and  operation  of  4th- order  tensors  becomes  a 
practical  necessity.  Several  familiar  examples  of  technical  4thr 
order  tensors  are:  elastic  compliance,  electrostriction,  2nd-order 
term  of  non-linear  thermal  expansion  and  environmental  swelling. 
There  remain  many  other  important  directional  dependent  physical 
properties  such  as  strength  and  certain  transport  phenomena  whose 
tensorial  characteristics  have  yet  to  be  adequately  verified 
experimentally.  In  order  to  facilitate  the  analysis,  measurement 
and  operation  of  such  physical  properties,  we  derive  the  invariants 
of  a  4th-order  tensor  and  suggest  a  geometric  interpretation. 
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Invariants  of  4th-0rder  Tensor 

Tensor  representation  of  physical  properties  which  are 
directionally  dependent  for  engineering  application  is  customarily 

in  Cartesian  coordinates.  For  this  reason  the  derivation  of 
invariants  of  a  4th-order  Cartesian  tensor  is  discussed.  To  avoid 

distinction  between  pseudo-invariants  and  invariants,  only  right- 
handed  Cartesian  coordinates  are  used.  The  invariants  derived 
are  for  rotation  about  the  axis.  Companion  invarinats  for 
other  axes  of  rotation  can  be  derived  through  simple  permutation 
of  the  indicies.  At  different  stages  of  the  derivation,  the 
symmetry  conditions  of  the  4th-order  tensor  are  used  to  simplify 
the  algebra.  Appropriate  modifications  which  are  cumbersome  but 
straightforward  must  be  made  for  skew-symmetric  tensors. 

We  consider  *  4th-order  symmetric,  Cartesian  tensor  and 

note  that  from  tensor  algebra  that  its  scalar  contractions  are 
invariant  to  orthogonal  coordinate  transformations.  We  shall  seek 
these  scalar  contractions  and  find  the  number  of  independent 
invariants. 

A  4th-order  tensor  can  be  contracted  through  combina¬ 

tions  of  the  substitution  tensor  6^^  and  the  permutation  tensor 

eijk 


which  are  defined  as: 
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100 

010 

001 


(1) 


and 


1+  1  for  even  permutation  of  ijk 
-  1  for  odd  permutation  of  ijk 
0  for  no  permutation 


The  2nd-order  tensor  components  of  the  4th-order  tensor  can  be 


obtained  by  contraction,  they  are: 


(2) 


5ik  Sijkl  = 


(3) 


Equations  (2)  and  (3)  are  the  only  distinct  components  because 
of  symmetry ,  i . e . , 


Sijkl  "  Sijlk 
Si jkl  *  Sklij 


(4) 


The  3rd-order  components  of  the  4th-order  tensor  are: 


emik  Sijkl  "  ®mjl 

emij  Sijkl  ”  ®mkl  (6) 

Again,  from  symmetry  conditions  Eq.  (4) ,  Eqs. (5  and  6)  are  the 
only  distinct  components.  Furthermore,  it  can  be  readily  shown 
from  the  properties  of  the  permutation  tensor  that  Q'j^  = 

The  2nd-order  components  of  Qmj^  then  takes  the  form: 


enjl  ^tijl 
Enkl  ®  mkl~ 


njl  Emik  Sijkl  Tnm 

(7) 

c  q  E  T '  =  0 

nkl  mij  bijkl  inm 

(8) 

0. 


First  and  Second  order  invariants  can  now  be  derived  for  those 
2nd-order  tensor  components  D,jl'  Tnm  ^ the  ^th-or^er 


tensor  S 


ijkl* 

The  first  order  invariants  are: 


4kl 

Dki  ’ 

6kl 

5id 

Si  jkl 

E  X1 

(9) 

6J1 

Dji  * 

6ik 

Si  jkl 

(10) 

6 

nm 

Tnm 

5 

nm 

Enjl 

emik 

Sijkl  5  X2 

(11) 

Expanding  Eq. (9)  we  can  express  the  first  order  invariants  of 

in  terms  of  its  components.  In  2-space,  i.e.,  for  i,j,k,... 
1,2,  they  are: 


I 


1 


Xi 


I 


2 


-  Sim  +  S2222  +  2S1122 

*  Sllll  +  S2222  +  2S1212 
='4S1122  +  4S1212 


(12) 

(13) 

(14) 


It  is  evident  from  Eq.  (12) , (13) ,  and  (14)  that  only  two  I's  are 
independent.  This  is  consistent  with  the  induction  that  there 
are  two  first  order  invariants  in  2-space  for  a  4th-order  tensor. 
\Te  arbitrarily  assign  I ^  and  I2  to  be  the  first  order  invariants 
to  be  consistent  with  underived  definitions  given  by  Hearmon^. 
Also,  to  be  consistent  with  the  definition  in  Ref.  [1],  Eqs. 

(11  and  14)  have  been  multiplied  by-2. 
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The  second  order  invariants  which  are  hitherto  not  known 
can  be  arrived  at  in  similar  manners.  By  taking,  what  is  in 
effect,  die  dot  products  of  the  cross  products  of  the  vector 
components  of  the  2nd-order  tensor  components  Dji'  Tnm j 

of  the  4th-order  tensor  the  following  2nd-order 

invariants  can  be  arrived  at. 


pq 

epik 

eqji 

T. 

13 

Tkl  " 

(al  ‘ 

-  2a2  + 

pq 

epik 

eqji 

D. 

13 

Dkl  = 

al  " 

a6 

pq 

epik 

eqji 

D!  . 
13 

Dkl  = 

a4 

a7 

pq 

epik 

eqji 

D.  . 
13 

Dkl  = 

a2  - 

a8 

pq 

epik 

eqji 

“u 

Tkl  = 

2  (a6 

-  “a) 

pq 

epik 

Eqjl 

Dij 

Tkl  = 

2  (°8 

*7) 

(16) 

(17) 

(18) 

(19) 

(20) 


^  can  be  expressed  in  terms  of  by  using  Eqs.(2,3  and  7) 

For  example: 

^1  -  ^pq  Epik  Eqjl  Eimg  Ejnh  Smngh  Eksu  Eltv  Sstuv  ^ 

=  6  e  e  <5  S  .  .  6  S  ,  , 

v2  pq  pik  qjl  mn  mni]  uv  uvkl 


For  conciseness,  the  remaining  of  \fj^  are  not  expanded.  The  ou 
expressed  in  terms  of  components  of  in  2-space  are: 


(16*) 
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“l  Si  i  j  j  S  k  k  i  1  bllll  +  4S1111SU22  +  2S1111S2222  +  4S1 122  +  4S1122S2222  4  S2222  (21) 


°2  '  Siijj  Sklki  SU11  4  2S1111S1122  +  2S1111S2222  4  2S1111S1212  4  4S1122S1212  (22) 

2 

f  2S1122S2222  4  2S2222S1212  4  S2222 

°3  *  Sijkl  Sikjl  *  Sllll  4  4 S 1 1 1 2  4  4S1122S1212  4  4SL22  4  2r5212  +  S2222  (23) 


“4  "  Sijij  Sklkl  "  Sllll  4  4S1111S1212  4  2SilllS,222  +  4S1212  4  4S1212S2222  +  S2222  (24) 


+  ^ini^nn  +  n  i®'  +  +  S_, 


a5  ■  Si jkl  Sijkl  *  Sllll  4  4sill2  4  2S1122  4  4S1212  4  4S1222  4  S2222 

®«  *  Siikl  S j jkl  *  Sllll  +  2S1111S1122  4  2 S 1 1 2 2  4  2S1112  +  4S2212S1112  +  2S2212 

4  2H1122S2222  4  S2222 

®7  "  Sikil  S j k j 1  *  Sllll  +  2S1111S1212  4  2S1112  +  4S1112S2212  +  2S1222 

4  2S12i2S2222  ♦  2  S 1 2 1 2  +  S2222 


aD  *  S  .  .  .  S 


8  '  “iijk  "ljlk  '  °1111  x  111  °1212  4  SllllSli''2  4  2S1112  +  2S1122S1212 


4S1112S2212  +  2S1 222  4  S1122S2222  4  S1212S2222  4  S22  22 
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It  can  be  shown  that  both  tl^  and  are  invariant  to  coordinate 
transformations.  Just  as  in  the  first  order  scalar  contractions, 
not  all  of  the  2nd-order  scalar  contractions,  4^  or  a^,  are 
independent.  Gauss-Jordan  reduction  revealed  that  the  system  of 
equations  Eq.  (21)  to  (28)  is  of  defect  three.  The  three  residual 
equations  are: 


“l  ‘  “2  ‘  2cx6  +  2a8  =  0 

a2  “  a4  +  2a7  “  2a8  =  0 

a3  "  a5  +  a6  +  a7  "  2a8  =  0 


(29) 

(30) 

(31) 


In  addition,  there  exist  three  relations  between  a's  and 
the  first  invariants  I1#I2  as  can  be  shown  from  Eqs.  (12) , (14) , (21) , 
(23),  (24)  and  (25). 

i 

X2  =  8(“5  ”  “3) 

(lx  +(1/2)I?)2  =  a4 

Equations  (21)  through  (34)  indicate  that  out  of  the  eight  a's 
only  two  a's  are  independent.  This  is  again  consistent  with  the 
mathematical  induction  that  there  are  two  2nd-order  invariants 
for  a  4th-order  tensor  in  2-space.  Any  two  a's  can  be  arbitrarily 
assigned  as  2nd-order  invariants.  We  chose  a3  and  a7  and  define 
the  following  combinations  as  2nd-order  invariants. 


(32) 

(33) 

(34) 


1 


I 

I 

I 


-  8 

«1  =  2a?  -  (i,  +  (1/2)  I2)2 

II2  -  8(0.3  -  c7)+  (13  +  I2)2 

We  will  show  in  the  following  section  that  these  particular 
combinations  lend  themselves  to  a  convenient  geometric  inter¬ 
pretation  akin  to  the  Mohr's  circle  for  2nd-order  tensors. 
Summarizing,  we  have  derived  1st  and  2nd-order  scalar  contractions 
for  a  4th-order  tensor.  We  also  show  that  there  are  two  independent 
lst-order  invariants  and  two  independent  2nd-order  invariants  for 
a  4th-order  tensor  in  2-space.  They  can  be  expressed  in  terms 
of  the  components  of  the  4th-order  tensor  components  as : 


I1  =  Sllll  +  S2222  +  2S1122 

(37) 

I2  *~4S1122  +  4S 

1212 

(38) 

11 1  =  (Sllll  "  S 

2222  j2  *  4(S1112  +  S2212  ) 

(39) 

XI2  "  (Sllll  +  S 

2222  "  2S1122  "  4S1212^  +  16^S1112 

S2212)2  (40) 

In  contracted  notation* , the  equivalent  expressions  are: 

I1  ■  S11  +  S22  + 

2®12 

(37') 

X2  ='4S12  +  S66 

(38') 

XI1  *  (Sll  -  S22 

f  +  (S16  +  S26)2 

(39') 

II2  =  (Sll  +  S22 

2S12  “  S66  j2  +  4  (S16  '  s2«)2 

(4'*) 

*The  contracted  notation  follows  the  common  practice  of  representing  a 
4-th  order  tensor  with  two  indexes  where  S^^z  Sn,  Si]  22  =  S12, 

2s1112  =  Si6 ,  2S2212  5  s26'  4S1212  E  s66*  For  more  detail  see  Ref.[l]. 

(35) 

(36) 
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Geometric  Interpretation  of  4th-0rder  Tensor  Invariants 


Because  of  the  large  number  of  components  involved,  the 
geometric  interpretation  of  tensors  of  order  higher  than  two  is 
a  formidable  task.  Even  in  the  plane  case,  the  geometric 
representation  for  a  2nd-order  tensor  is  equivalent  to  represent 
in  two  dimensions  a  system  of  three  quadratic  trigonometric 
functions.  Mohr's  circle  provides  such  a  representation.  It 
contributed  to  better  visualization  of  a  2nd-order  tensor  and 
provided  greater  operational  convenience  for  engineering  problems 
in  the  analysis  of  stress  and  strain  as  well  as  for  problems  in 
rigid  body  dynamics.  Similarly,  the  representation  of  a  4th-order 
tensor  is  equivalent  to  representing  a  system  of  4th-order 
trigonometric  functions.  A  geometric  representation  similar  to 
that  of  the  Mohr's  circle  has  been  suggested  by  P.  Mast  [2], 
Utilizing  the  invariants  derived  in  the  previous  section  Eqs.  (37), 
(38) ,  (39)  and  (40),  this  can  be  generalized  to  geometrically 
represent  the  transformation  of  a  4th-order  tensor  from  any 
arbitrary  non-principal  direction. 

From  the  definition  of  tensors,  for  rotation  about  the  x^  axis, 
the  components  of  a  4th-order  tensor  at  an  arbitrary 

orientation  0  with  respect  to  a  material  coordinate*  can  be 
transformed  to  the  components  at  another  orientation  0' 

through  the  relation  (1] : 


♦Without  loss  of  generality,  we  may  consider  the  material 
coordinate  coincides  with  the  principal  direction. 
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,,  i  r 


;;  2 


s  ' 

1122 


S'  r 
2222 


m  2  2  4  4 

m  n  m  tn 


mn  -ns  r.  m2n2  m^n-mn^ 


■)  ■> 
-m"n 


Si 


3n  2  (m  n-mn3)  m4-3m2n2  2mn3  3m2n2-n4  mVmn3  i  .IS 


4  _  2  2 

n  2m  n 


-2m3n 


2  2  i 

m  n  | 


2212 ^  !  -2mnJ  2 (mn J-m Jn )  3m2n2-n4  2m3n  m4-3m2n2  mn3-m3n 


4Si2l2J  ^4rr,2n2  4(mn3-mV  4m2n2  4  (in3r.  -mn3)  (m2-n2) 

where  m  and  n  are  directional  cosine  and  sine. 


33)  ._2_2 


2212 


b"'21 


The  graphical  representation  of  this  4th-ord>r  tensor  transformation 
can  be  observed  by  rewriting  Eq .  41  in  multiple  c,ngle  representation 
as  described  by  Tsai  and  Pagano  [3]  in  the  following  form: 


where 
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4  *  e  -  e1 

u!  -  V8  (ai,  +  i2 ) 

U2  =  1/2  |  S2222  "  Sllll) 

U3  =  1/8  ^  4S1212  +  2S1122  "  Sllll  "  S2222) 

U4  -  1/8  ( I,  - 

U5  =  1/8  (*2  +  *l) 

U6  =  (S1112  +  S2212^ 

U?  =  1/2  ( S1112  -  S2212j 

The  geometric  interpretation  of  this  tensor  transformation 
can  be  recognized  after  some  rearrangements.  Take  the  first  of 
Eq.  (42)  for  example, 

“  U1  ”  U2  cos  2<^  +  U6  s*n  2<^  “  U3  COS  4<^  +  U7  s*n  4<^ 
Consider  the  terms  containing  2<f>  and  define: 

U'  =  U2  cos  2$  -  Ug  sin  2$ 

If  U2  and  Ug  are  the  sides  of  a  right-angled  triangle  with  a 
hypothenuse  R^,  as  shown  in  Fig.  1,  then 

U2  =  Rj^  cos  26 
Ug  =  R1  sin  26 

R1  =  U2  cos  20  +  ug  s^-n  20 


(43) 


(44) 


(45) 


(46) 


Substitute  Eq.  (46)  into  (45),  we  obtain 
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U2  =  *U2  COS  26  +  U6  Sin  20)cos(20  +  2<f>)  <47) 

=  cos  26' 

Similarly,  we  can  express  terms  containing  4<J>  as: 

cos  4<J>  -  sin  44>  =  R2  cos  40*  (48) 


where 

R2  =  U3  cos  44>  +  U7  sin  4<j> 

Making  use  of  Eq.  (47)  and  (48),  Eq.  (44)  simplifies  into  the  form: 


Silll  =  ui  “  Ri  008  201  “  r2  cos  401 


(49) 


As  shown  in  Fig.  2,  Eq.  (49)  can  be  represented  by  the  horizontal 
projected  distance  between  the  radius  vectors  of  two  circles 
separated  by  distance  between  centers  and  with  radii  and  R2 
respectively.  The  radius  vector  for  R^  rotates  at  20'  and  that 
for  R2  rotates  at  40'.  Both  the  original  at  0  and  the 

transformed  at  0'  are  shown  in  Fig.  2.  It  can  also  be 

easily  shown  that  R^  and  R2  are  the  roots  of  the  2nd  invariants 
[Eq.  (39),  (40)]. 


Carrying  out  similar  rearrangements,  Eq.  (42)  can  be  written  as: 


• 

-  13  - 

C  I 

bllll 

U1 

-R1 

0 

-R2 

0 

1 

fa1122 

"U4 

0 

0 

R2 

0 

cos  20' 

b1112 

0 

0 

1/2RX 

0 

R2 

sin  20  ' 

C  1 

b2222 

= 

U1 

R1 

0 

”R2 

0 

cos  40  1 

b2212 

0 

0 

i/2Rl 

0 

CM 

« 

1 

sin  40' 

C  1 

S1212 

_U5 

0 

0 

-R2 

°J 

0 

L  -1 

The  components  of  the  4th-order  tensor  at  orientation  0  (which  is 
equivalent  to  $  *  0  in  Eq.  51)  can  be  represented  geometrically 
in  Fig.  3.  In  the  construction  of  this  representation,  it  is 
assumed  that: 


U2  >  0  or  S2222  >  Sim 

u3  >  °  or  (4S1212  +  2sH22)  >  (Sllll  +  S2222^  (52) 

U4  >  0  °r  4s1212  >  (  Sllll  +  S2222  +  6S1122j 


In  addition,  just  as  in  the  Mohr’s  circle,  certain  sign 
conventions  must  be  followed.  A  positive  rotation  <j>  (counter¬ 
clockwise)  in  the  material  plane  is  equivalent  to  a  positive 
rotation  of  2$  (counter-clockwise)  of  radius  vector  and  a 
positive  rotation  4$  (also  counter-clockwise)  of  the  radius 
vector  R2*  The  signs  of  the  components  are  determined  by 

the  direction  of  the  measurements.  All  measurements  from  left 
to  right  and  those  from  bottom  to  top  are  positive.  In  order 
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to  avoid  confusion,  the  directions  of  measurements  as  indicated 
by  the  arrows  in  Fig.  3  must  be  followed.  For  example, 
is  obtained  by  the  vertical  distance  from  the  tail  of  the 
diameter  vector  of  circle  R2  to  the  head  of  diameter  vector  of 
circle(l/2)R^.  Note  that  when  0  >  v/2,  the  measurement  direction 
for  S1^12  is  reversed  and  takes  on  negative  values. 

The  above  geometric  representation  can  be  constructed  from 
any  given  set  of  Si;jkl.  The  orientation  of  the  principal 
direction  need  not  be  known.  The  invariants  1^  I2,  rI2 

can  be  computed  from  Eqs.  (37,38,39  and  40).  From  Eqs.(43  and  50), 
U  ij4,  u5,  R1  and  R2  can  be  computed.  These  values  can  be  used 
to  locate  the  centers  of  the  circles,  reference  lines  A  and  B 
and  the  circles  , (1/2) R1# and  R2  as  shown  in  Fig.  3.  Either 
S  or  S1010  can  then  be  used  to  determine  the  orientation  of 
the  radius  vector  R2<  Thus  the  principal  direction  0  can  be 
directly  determined,  which  in  turn  determines  the  radius  vector  R^. 
The  components  of  the  4th-order  tensor  S|jk^  at  orientation  <J> 
from  Si;.kl  can  be  determined  by  rotating  the  R1  by  2<J>  and  R2  by  4$ 
and  the  appropriate  distances  measured  according  to  the  suggested 
sign  conventions. 

The  representation  presented  herein  is  but  one  of  many 
available.  It  serves,  primarily,  to  illustrate  the  essential 
features  which  exist  in  the  geometrical  representation  of  a 
4th-order  tensor.  Detail  variations  will  depend  on  personal 
preferences  and  specific  application.  For  example,  if  the 
transformation  of  Si;jkl  in  contracted  notation  is  desired,  the 
third,  fifth  and  sixth  equations  of  Eq.  (51)  must  be  multiplied  by 
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2,  2  and  4  respectively.  It  follows  that  in  the  geometric 
representation,  the  circle  (1/2)^  can  be  removed  and  concentric 
circles  2R2  and  4R2  added  to  the  cricle  Rj.  In  addition,  the 
location  of  reference  lines  A  and  B  as  well  as  the  relative 
positions  of  the  circles  are  arbitrary.  If  they  are  changed 
from  the  configuration  presented  here,  appropriate  sign  convention 
has  to  be  re-established. 

It  is  also  worthwhile  to  note  the  analogy  between  the  Mohr's 
circle  transformation  for  2nd-order  tensors  and  the  representation 
for  4th-order  tensors.  In  Mohr's  circle,  the  first  order 
invariant  determines  the  location  of  the  center  of  the  circle , 
and  the  2nd-order  invariant  determines  the  magnitude  of  the  circle. 
In  the  4th-order  representation,  the  first  order  invariants 
determine  the  location  of  the  centers  of  the  circles  and  the 
reference  lines  A  and  B  while  the  2nd-order  invariants  determine 
the  magnitudes  of  the  circles.  For  a  rotation  tj)  in  the  physical 
plane,  the  second  order  transformation  is  represented  by  a 
rotation  2<J>  in  the  Mohr's  circle  while  the  4th-order  transformation 
xs  represented  by  rotations  of  2$  and  44>  of  the  circles. 

In  practical  application,  this  geometric  representation  is 
used  to  graphically  transform  the  compliance  matrix  S^..^  of 
fiberglass  reinforced  composites.  It  was  found  that  upon  the 
computation  ->f  the  invariants  and  the  construction  of  the  circles, 
any  jkl  can  be  rapidly  determined  by  laying  out  two  vectors 
r  and  R_  at  the  desired  orientations  20  and  40.  From  constructions 

X  “ 

on  8"x  10 "graph  papers*, the  accuracy  attained  is  within  1%  of 
the  computer  computed  results  using  Eq.  (41)  . 

♦Graphical  construction  can  be  simplified  by  normalizing  Eq.(51)  by 
Rx  pr  R2  or  U]_, 
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In  cases  where  principal  direction  is  to  be  found,  it  can  be 
determined  graphically  by  the  intercept  of  S{212  to  c^rc^e 

,  and  this  is  found  to  be  much  faster  than  analytically  setting 

S|ii2  or  S 2212  to  zero  *-n  and  so^-ve  f°r 

Applications 

We  have  derived  the  first  and  second  order  invariants  for 
4th-order  tensors  and  have  shown  that  in  2-space  there  exists 
two  independent  first  order  invariants  and  two  independent 
2nd-order  invariants.  By  expressing  the  4th-order  tensor 
transformation  in  multiple  angle  representation  and  making  use 
of  the  invariants,  a  geometric  representation  analogous  to  the 
Mohr's  circle  can  be  constructed. 

The  invariants  and  graphical  representation  presented  here 
are  applicable  to  4th-order  tensors  encountered  in  the  physics 
and  mechanics  of  solids  and  fluids.  They  are  also  applicable  for 
certain  engineering  constants  for  composite  materials  which 
transform  in  the  form  of  Eq.  (1).  Thus,  similar  invariants  and 
geometric  representations  can  be  derived  for  the  A,  B,  and  D 
matrix  for  laminate  plates  [4].  Familiar  applications  of  the 
invariants  of  2nd-~order  tensors  and  Mohr's  circle  suggest  some 
natural  applications  of  their  counterparts  for  4th-order  tensors. 
For  example,  the  principal  values  and  principal  directions  of  the 
stiffness  tensor  for  crystals  and  composite  materials  can  be 
conveniently  obtained  from  the  invariants  and  the  graphical 
representation.  In  view  of  the  contribution  of  stress  invariants 
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-o  the  formulation  of  yield  criterion  for  isotropic  materials,  the 
4th -order  tensor  invariants  may  be  useful  in  the  characterization 
of  the  flow  and  fracture  of  anisotropic  solids.  Finally,  the 
newly  derived  second  order  invariants  may  be  considered  as 
additional  intrinsic  material  properties  and  may  be  used  to 
simplify  optimization  of  the  physical  properties  of  laminated 
composites  as  suggested  in  Reif.  [3]. 
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Appendix 

Invariants  of  4th-order  tensor  in  3-space : 

The  derivation  cf  the  invariants  of  a  4th-order  tensor  in 
3-space  for  a  rotation  of  coordinate  around  the  X^-axis  is 
similar  to  that  in  2-space.  Because  X^  =  x^  in  the  transformation, 
the  components  of  which  contain  the  3  component  transform 

as  tensors  of  orders  less  than  four;  i.e.,  the  components 


S1123, 

S2213, 

S1113, 

S2312, 

S2223 

S1312 

Transforms 

as 

3rd-order 

tensor 

(Al) 

'1313' 

S2323' 

S1323  j 

|  Transform 

as 

2nd-order 

tensor 

(A2) 

1133' 

S2233' 

S3312  j 

3323' 

S3313 

Transform 

as 

lst-order  tensor 

(A3) 

Thus,  in  addition  to  the  invariants  derived  (Eq.  37  to  40),  there 
exist  invariants  associated  with  the  tensor  components  in  Eq. (Al  to  A3). 
For  the  components  in  Eq. (Al) ,  it  can  be  easily  shown  that  there 
exists  no  lst-order  invariants.  The  2nd-order  invariants  are: 


”3  '  1 

^1123  ~  S2223  +  2S1312)2  + 

(S1113  '  S^213 

2S2Zu)2 

(A4) 

“4  "  1 

(3'j1123  +  S2223  ’  2S1312  )  2 

+  | S1113  +  3S2213 

2S2312j 

(A5) 

II5  "  1 

(3S1113  +  S2213  +  2S2312)  + 

(S1123  +  3S2223 

+  2S1312)2 

(A6) 

“6  ’  ( 

®1113  “  S2213  +  2S2312)2  + 

(S1123  -  S2223 

2S1312)2 

(A7) 

For  the  components  in  Eq.  (A2) ,  the  first  and  second  order  invariants 
can  be  obtained  from  well  known  results  of  2nd-order  tensor 
transformation.  The  first  order  invariants  are 
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r3  * 

4S2323  +  4S1313 

(A8) 

S2233  +  S1133 

(A9 ) 

The  second  order 

invariants  are 

ii 

!"• 

H 

M 

(S1133  "  ^2233  j 

!  +  (JS3312)2 

(A10) 

H 

M 

CO 

II 

( S2323  "  S1313  ) 

+  (4!i2312  \ 

(All) 

Finally,  for  the  components  in  Eq. (A3) ,  no  first  order  invariants 
exist  and  the  second  order  invariant  is  the  magnitude  of  the 


vector 


(A12) 


Recapitulating,  for  a  4th-order  tensor  in  3-space  in  a 

rotation  around  the  X-axis,  the  lst-order  invariants  are:  I^Eq.37), 


I2(Eq.38),  I3(Eq.A8),  I4(Eq.A9), 

and  the  2nd  order  invariants  are  II^Eq.38),  II2(Eq.39),  IT3(Eq.A4), 
II4(Eq.A5),  II5(Eq.A6),  II6(Eq.A7),  II7(Eq.A10),  II8(Eq.All), 

IIg  (Eq . A12) . 


The  geometric  representations  of  the  transformation  of  the 
components  in  Eq.  (A2)  and  Eq. (A3)  are  that  of  the  conventional 
Mohr ' r  circle  and  that  of  the  vector  circle  respectively.  The 
geometric  representation  of  the  components  in  Eq. (Al)  can  be 
seen  by  expressing  the  transformation  in  multiple  angles  using 
the  relationships 


sin3  6  =  1/4(3  sin  0  -  sin 
2 

cos  0  sin  0  =  1/4  (cos  0  - 
2 

cos  0  sin  0  =  1/4 (sin  0  + 

3 

cos  0  =  1/4(3  cos  0  +  cos 


30) 

cos 

sin 

30) 


30) 

30) 


(A13) 
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q  t 

b1123 


R.  0  R-  0 

4  3 


COS  0 


q  l 

b1113 


Rc  0 

3 


sin  0 


q  i 

2223 


r5  0  -r3  0 


cos  30 


(A1 4  ) 


q  I 

°2213 


R  a  0  -R, 

4  3 


sin  30 


q  ' 

2312 


■R6  0  -R3 


q  < 

b1312 


-Rg  0  R3  0 


wbare 


R3  "  l/2(lI3)V2.  R4  -  J/2(lI4)1/2,  R5  -  l/2(ll5)1/2. 


R6  =  1/4 


(K«): 


The  geometric  representation  of  equation  (A14)  is  similar  to  the 
one  constructed  for  Eq. (51)  and  is  shown  in  Fig.  4.  In  this 
construction,  II 3,  II^,  II5'  IJ6  are  assumec*  to  positive.  It 
can  be  noted  that  in  the  principal  direction  (0  =  0), 

S1113  =  S2213  =  S2312  =  0  indicating  that  the  material  is  symmetric 
to  the  2-3  plane,  i.e.,a  mono-clinic  material  with  as  the 
principal  axis.  It  is  important  to  note  that  since  the  transformation 
of  in  2-space  (Eq.  41)  is  not  coupled  to  the  transformations 

of  the  components  in  3-space  (Eqs.Al,  A2  and  A3),  the 

invariants  Eqs.  (37-40)  and  geometric  representation  (Fig.  3) 
remain  the  same  for  3-space.  Previous  remarks  on  the  applications 
of  the  invariants  and  geometric  representation  for  2-space  are 
also  applicable  for  their  counterparts  in  3-space. 
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FIGURE  2 


FIGURE  4.  TRANSFORMATION  OF  THIRD  ORDER  TENSOR 
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